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Abstract 

Crook’s Fluctuation Theorem (CFT) and Jarzynski equality (JE) are effective tools for obtaining 
free energy difference AF(Xa_ —>• A^, To) through a set of finite-time protocol driven non-equilibrium 
transitions between two equilibrium states A and B (parameterized by the time-varying protocol 
A(t)) at the same temperature Tq. Using a new work function AWg, we generalize CFT to transi¬ 
tions between two non-equilibrium steady states (NESSs) created by a thermal gradient and show 
that it is possible, using the same set of finite time transitions between these two NESSs, to ob¬ 
tain AF{Xa —>■ Xb,Tq) for different values of Tq, thus completely eliminating the need to make 
new samples for each new Tq. The generalized form of JE arises naturally as the average of the 
exponentiated AWg- The results are demonstrated on two test cases: (i) a single particle quartic 
oscillator having a known closed form AF, and (ii) a 1-D chain. Both systems are sampled from 
the canonical distribution at an arbitrary T' with A = A^, then subjecting it to a temperature 
gradient between its ends, and after steady state is reached, effecting the protocol change Xa ^ Xb 
in time r, following which AWg is computed. The reverse path likewise initiates in equilibrium 
at T' with X = Xb and the protocol is time-reversed leading to A = Ayi and the reverse AWg- 
Our method is found to be more efficient than either JE or CFT when free-energy differences at 
multiple To’s are required for the same system. 
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Consider a thermo-mechanical system whose equilibrium state is dehned by its tempera¬ 
ture To and an external protocol A fixed at (for examfoe the position of confining potential 
|l|, the position of the last molecule of a protein chain 2| etc.). A large class of problems in 
biological and chemical physics (such as transition between conformations of proteins, fold¬ 
ing and unfolding of proteins, enzyme-ligand binding, hydration etc.) concerns the change in 
free energy, AT(Aa — \b,To), of this system as its configurational space evolves under \(t) 
in a finite time r corresponding to the final value X = Xb and the system eventually relaxes 
to a new equilibrium at the same temperature Tq. Several methods have been proposed 
for computing AF{Xa —t A^Tq) - thermodynamic integration j^, umbrella sampling jd], 
steered molecular dynamics [5 1, and nonequilibrium work relations {gIJhI. 

The development of Jarzysnki’s equality (JE) j^8] and Crooks’ fluctuation theorem 
CFT) 9|, [lO, has dramatically improved our ability to calculate free-energy differences 
13- 3 of real systems js, 17| through finite-time irreversible processes between two equi¬ 


librium states at the same temperature, Tq. Nevertheless, the task remains daunting be¬ 
cause of the requirement of extensive sampling of the configurational space. In addition, 
AF{Xa —^ As, To) thus computed is valid only for the particular temperature at which the 
samplings are performed and if AF{Xa Xb,T' ^ Tq), is needed, the re-sampling of the 
entire data set is necessary at T'. 

In this work we generalize CFT and JE by proposing a new fluctuation theorem that 
enables us to calculate AF{Xa As,To) with good accuracy for a range of Tq values using 
a single set of sampling data, thereby completely eliminating the need to make new samples 
for each new Tq. The proposed fluctuation theorem utilizes the transition between two 
nonequilibrium states, and Tq features in the equation as a scaling parameter. 

Let us now look at details of the problem. For a system in canonical equilibrium, the 
Helmholtz free energy is: 


F{TQ,X) = -kBTQlog(^J exp[-/3oT(r,A)]dr^ , (1) 

where /c^ is the Boltzmann constant and (3q = (/csTq)”^. The system’s energy T(r, A) de¬ 
pends upon the microstate T and varies parametrically over time according to E{r(t), X(t)) = 
+ ^{.xi,X2, ■ ■ ■ ,xn, X{t)), where pi and Xi are the momentum and the position of 
the particle. In CFT, the system is initially in equilibrium state A with A = Ay^. At time 
t = 0, X starts to evolve until t = r, and stays fixed at its new value Xb- During this period 
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work XdE/dXdt is performed. The superscript F denotes the forward transition 

A ^ B. Over time, the system relaxes to a new equilibrium state B with X = Xb- Being 
irreversible, the work depends upon the initial microstate r(0) of the system (and its 
surroundings), and therefore, an exhaustive sampling of the initial microstates provides the 
probability density of forward work, P{W^ = w). Now consider the same system evolving 
in a reverse manner. The system begins at equilibrium state B where X = Xb, and over 
0 < f < r, A traces itself back from Ab —?■ A^. Eventually the system reaches the equilibrium 
state A. Repeated sampling of this reverse transition provides P{W^ = —w) for the reverse 
work. GET relates the ratio of these two densities with AF{Xa Xb,To): 

P{W^ = Wo) 


= exp[-/3o(wo - AF(To, A))] 


( 2 ) 


P{W^ = -Wo) 

The validity of [2] requires the dynamics to be ergodically consistent i.e. if a microstate has 
a nonzero probability in equilibrium state A, it evolves to a microstate that has a nonzero 
probability in equilibrium state B. Integrating ([2]) gives JE |^. However, since To is implicit 
in the sampling dynamics, the probability densities obtained cannot be used to calculate 
AF(r', A) if T' 7 ^ Tq. In order to employ a single set of sampling data for calculating 
AF{T', X) corresponding to a range of temperature T', the dependence of sampling data on 
To must be removed. We set out to do this by looking at the work and heat distributions 
during the transition between two nonequilibrium steady states. 

Rather than beginning at equilibrium, we begin at a nonequilibrium steady state SSi 
obtained by imposing a temperature difference (T^ — Tc) at the two ends of the conductor, 
where Tjj and Tc are the temperatures of the hot and cold ends. This steady state originated 
from some primordial arbitrary equilibrium state A characterized by A^ and Tq by employing 
suitable temperature constraints. For all practical purposes, the system reaches steady-state 
when the relevant time-averaged macroscopic observables become stationary. Th, Tc and Tq 
are related to each other through Th = Tq-I-ATb and Tc = To — ATc- Thus, depending upon 
ATh and ATc, both not necessarily being equal, one can think of starting from arbitrarily 
different canonical equilibrium states. Note that this allows us to choose any arbitrary Tq. 

After SSi is achieved, at f = 0, A starts to evolve from Aa until time t = t when 
X = Xb and work is performed. This external work does not result in any phase-space 
compression. Given sufficient time, the system reaches a new steady state S'S' 2 - Upon 
removing the temperature constraints, the system eventually reaches the equilibrium state 
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B, defined by As and Tq. The reverse transition can likewise be accomplished under the 
time-reversed protocol. Such transition between steady-states has been studied before in 
a different context 19|] . The underlying principle governing our approach is the relaxation 


of a nonequilibrium state to an equilibrium state 


20| . This relaxation is governed by the 


constraints imposed on the system, and thus one can obtain a multitude of equilibrium states 
from a single nonequilibrium state by judiciously choosing the constraints and boundary 
conditions. 

In state A (state B), the system follows the canonical distribution parameterized by 
(by Ab): 

/e,,A(r) = ^ exp [-(3oE{T, A^)] (3) 

The density function of the nonequilibrium state and the Jacobian are given by Liouville’s 


equation 


M 


/55,l[r(t)] = /eg,Al [r(0)]exp 


J dt’{Anit') + Ac{t')) 


dT{t) = (jr(O)exp 


/ dt' (Ab^(T) -|- Ac{t')) 


to 


( 4 ) 


where A = 


dV/dV denotes the phase-space compression factor, with H {C) denoting the 
hot (cold) region. The intermediate region does not contribute to A (owing to Hamilton’s 
equation of motion). Importantly, the normalizing constant corresponding to SSi is the 
same as the partition function for A. The phase-space compression factors are related to 


the heat flow 


22 


25| from the thermostats through: 


(Qjr)t — kBTH{AH)tt , {Qc)t — kBTc{Ac)tt- 


( 5 ) 


For sake of compactness, we will drop t from the density functions and cumulative heat 
flows later. Next, we bring the generalized dimensionless time-integrated work function, 
AWcit) 26| into picture, which can relate two microstates (r(0) and r(t)), neither of them 
necessarily in equilibrium: 


exp {AWcit)) 


/i(r( 0 ))dr( 0 )ZAo 
um)dT{t)z,, ’ 


( 6 ) 


The initial microstate r(0) evolves to r(f) in time t. /i(r(0)) (or f 2 (T(t))) is the probability 
density of r(0) (or T(t)) corresponding to an associated equilibrium state 1 (or 2). We 
conjecture that such an association is possible after the system undergoing non-equilibrium 
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transition loses its memory. Z\^ denotes the partition function at A*. Now we bring the 
superscripts F (for the forward transition A ^ SSi ^ SS 2 ^ B) and R (for the reverse 
transition B —)■ SS 2 —)■ SSi —?■ A). The forward transition takes r(0) — ?■ r(t), while the 
reverse transition takes r*(0) —)■ T*{t), where r*(0) is related to r(t) through time-reversal 
mapping. The generalized work function during A ^ SSi is (see Section-I of Appendix): 

1 /■* 


^^G,A^SSi 


Tr 


0 Jo L 


ATh ATc 

VH rr '^C 


dt' 


(7) 


Th Tc 

Proceeding analogously (see Section-Ill of Appendix), the generalized work function during 
A SS 2 is: 


^^G,A^SS2 


tI 


t-\-r 


/^F 


Qh- 




dt' 


Th Tc 

Therefore, the work function during SSi -A- SS 2 can be obtained by subtracting [7| from [SJ 


( 8 ) 


+ TfT 
In 


ATh ^F _ ^G , 
T 


H 


-Qkr - 


c 


(9) 


where the heat flows are for the time duration r over which A changes. In a similar manner, 
we can compute the work function during the reverse transition SS 2 —)■ SSp 




R 


rri ^H,T 

J- r 


H 




'G 


( 10 ) 


Now we make the important assumption of the ergodic consistency being valid during 
the transition SSi -A- SS 2 , and therefore, using [6] we can write: 

/eg,A [^SS,1 (0)] dTss,i{^)^XA 


exp — 


( 11 ) 


feq,B [rSS,2('r)] drsS, 2 {T)Z\g 

The subscripts SSi emphasize that the points are on trajectories whose evolution is described 
by equations of motion that take the ensemble of states from S'S'i at 0 to SS 2 at r. Because 
of the deterministic nature of the dynamics, AITg = —AWq For simplicity, 

we now drop all subscripts except G. The probability densities of the forward and reverse 
work functions therefore can be related as (see Section-II of Appendix): 


P [AW§ = -k]= e-^^P [AW^ = k] 


-‘Xb 


( 12 ) 


A rearrangement results in the proposed fluctuation relation-. 


P 


PoW" + A 




■H 


G 


= k 


P 


+ flo 


ATh^^ ATc 


’-H 


Qfl- 


'-c 


Q 


= -k 


^ g[fc-/3oAF] 


(13) 
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which is the main result of this paper (henceforth, referred to as GCFT). Since the samplings 
have been performed at Th and Tc, the effect of Tq is inherently absent in them, and 
is simply a scaling parameter. Depending upon the temperature at which AF is to be 
calculated, we can compute the forward and reverse densities of the work function simply 
by substituting the desired value of Tq. The generalized JE may be obtained by averaging: 


-yow^-fio 


ATh ^ ATc ^ 

i ( 


T 


H 


'-C 


^ g[-/3oAF]_ 


(14) 


A second law type inequality can be recovered by applying the Jensen’s inequality to ffTTD : 

( 16 ) 


AT AT 


Th Tc 

It must be noted that the above equations are not exact relationships, and hold true only 
for large r. Taking r large enough, while fixing the time required to reach the steady state, 
ensures that the contributions arising from phase-space compressions become negligible. We 
test the effectiveness of (1131) on a 1-D 0"^ chain of N particles. Its energy function is: 


N 


E = Y. + YU {Xi, Xi+i) + ^ E {Xi) . 


i=l 

\2 


N 

E 

2=1 


(16) 


i=i ‘^rui 

Here U{xi,Xi^i) = 0.5/ci (|xj+i — Xi\ — dY represents the quadratic nearest neighbour in¬ 
terparticle interaction, while V{xi) = 0.25/c2 (a^i — a(i,o)^ represents the quartic tethering 
potential with Xi^ being the equilibrium position of the particle. We have kept ki = 1.0 
and rui = 1.0. /c 2 plays the role of A: 


k2 = 0.25 Vt < 0 

= 0.25 (1 + lOt/r) VO < f < r 
= 2.75 Vf > r 


(17) 


Test Case 1: The first test case involves a single particle system (subscript 1 dropped) 
having a known analytical solution for AF\ 


AF = -fceTolog 




J e\--Po^^ e\--hop'll‘A dxdp J (18) 

= kBTo/4log(ll) 

We compare this known AF with our results. We subject the single quartic oscillator to a 
position-dependent temperature field. 


T(x) = 1-1- 0.1 tanh(a:). 


(19) 
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to bring it away from equilibrium. Temperature is controlled by Hoover-Holian thermostat 
j^ . The system is simulated for 100,000 time steps (each time step = 0.001) under this 
temperature held through which it reaches SSi. k 2 is changed over the next r = 10,000 time 
steps according to equation flTT)) . AVTg in this case is: 


AITg = (3oW + ^0 / Qdt' - / l3{x)Qdt\ 


( 20 ) 


where, (3{x) = l/ksT^x) and Q = —r}T{x) — 3p‘^^T{x). t] and ^ are the Hoover-Holian 
thermostat variables. W = 'H(r) — 'H(O), where 'H(t) = p^/2 + V + J{r]T{x) + 3p^^T{x))dt\ 
Probability densities of generalized work are constructed using 60,000 random initial points. 
Figure [T] shows AF due to the evolution of k 2 as a function of temperature: GCFT is able 
to reproduce the theoretical results accurately for a range of temperatures without the need 
to resample at every new Tq. 

Test Case 2: We now consider a larger system {N = 25). The system is initialized 
with Xi = Xifl = i and random particle velocities. The equations are integrated using 
classic Runge-Kutta algorithm with an incremental time step of 0.01. Post initialization, a 
temperature gradient is imposed on the systeni by keeping the two end particles at Th and 
Tc using two Nose-Hoover (NH) thermostats 28j]. Subsequently, after 1 million timesteps 
(steady state is assumed to have reached), k 2 evolves in r = 100,000 time steps. The 
cumulative heat flow from the hot thermostat is Qh = — jlTHPudf (likewise for the cold), 
where rinivc) is the hot (cold) NH variable. The work done due to the change in tethering 
potential during time r is hF = (r) — (0) where 


'H{t) = + ^ + St ricPc'it' + / riHPudf 


( 21 ) 


Here pn ipc) denotes the hot (cold) particle’s momentum, and $ = + ■ AWq 

and AWq are computed using 5,000 trajectories each. Figure [2] shows probability densities 
of the forward and reverse generalized work functions P{AWq) and P{AWq) at Tq = 0.29. 
Two pairs of {Th, Tq) - red for (0.27,0.23) and blue for (0.30, 0.20) - are chosen. The points 
of intersection of the forward-reverse pair gives PqAF flT^ which should be independent of 
{Th, Tq) for the same Tq as evident from the hgure. Importantly, these same 10000 samples 
can be used to compute AF at any Tq. Table ([T]) shows seven such Tq values, computed 
using both sets of {Th,Tc). Not only is AF at a given Tq independent of {Th,Tc) as it 
should be, it is clear that Tq does not even need to be within the range of {Th,Tc) for the 
method to work. 
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FIG. 1. Comparison of AF obtained using theoretical and proposed approaches in test case 1. 
Notice, that the proposed approach provides a good approximation to the theoretical results. 


Finally, Table (Plists AF computed using JE and CFT at the seven different tempera¬ 
tures. While GCFT is able to identify AF as accurately as CFT and JE, it does so with only 
one set of samples. CFT and JE on the other hand would require a new set of samples for 
each To, thereby imposing a severe computational or experimental burden on the analyst. 
We must, however, point out that the transition needs to be carried out slowly, as our efforts 
to calculate AF using r = 100 steps did not yield any fruitful result. 

To summarize in this work, generalized versions of CFT and JE have been presented. 
The proposed extensions present a suitable method through which equilibrium free energy 
differences can be extracted from the information embedded within the non-equilibrium 
steady states. The augmented equations bear remarkable similarity with those of CFT and 
JE with additional contributions arising due to heat flowing from the reservoirs. CCFT has 











Generalized Work (AW_) 

(j 

FIG. 2. Forward and reverse probabilities of generalized work function at Tq = 0.29 calculated 
using 5,000 forward and reverse trajectories (red) {Th,Tc) = (0.27,0.23) and (green) {Th,Tc) = 
(0.30,0.20). The forward and reverse probabilities at approximately the same value of Wq- AF 
calculated compares well with that out JE and CFT. Results obtained using the same dataset for 
other To values are similar, and agree well with CFT. 

been tested using two different cases, with each of them suggesting that GCFT is a suitable 
alternative to CFT and JE when evaluating AF at multiple temperatures. 


I. APPENDIX 

A. Section-I 

In this section we will derive equation (7) of the manuscript. Let us look at a system 
initially in canonical equilibrium (state A, temperature Tq), whose distribution function 
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TABLE I. Comparison of free energy differences using JE, GET and GCET for seven different 
values of Tq. GCFT results are for two different steady-states: Th = 0.27, Tc = 0.23 and Th = 
0.30, Tc = 0.20. Notice that the AT obtained using GGFT matches closely with those from JE 
and GET. It is interesting to note that the case of Th = 0.27, Tc = 0.23 is able to approximate 
the equilibrium free energy differences even for the states as far as To = 0.21 and Tq = 0.29. The 
results indicate that one can use a single set of data obtained during a transition between two 
NESS and employ GCFT to calculate free energy differences for a range of temperature. 


To 

JE 

GET 

GGFT 

{Th,Tc) = (0.30,0.20) 

GGFT 

{Th,Tc) = (0.27,0.23) 

0.21 

1.39 

1.35 

1.40 

1.39 

0.22 

1.48 

1.50 

1.49 

1.46 

0.24 

1.64 

1.60 

1.67 

1.62 

0.25 

1.73 

1.70 

1.75 

1.69 

0.26 

1.81 

1.79 

1.83 

1.78 

0.28 

1.98 

1.97 

2.02 

1.97 

0.29 

2.07 

2.00 

2.10 

2.07 


given by: 

/e,,^(r(0)) = (22) 

for a microstate r(0) of A. Za is the partition function and (3o = {kBTQ)~^. On this 
system, we apply a temperature gradient by keeping the two ends at temperatures Th 
and Tc- Because of the dynamical nature the system evolves to a new microstate r(t) in 
time t under the influence of the thermal gradient. Assurmng deterministic dynamics, the 
system evolves according to Liouville’s continuity equation 29|], and the evolved distribution 
function becomes: 


= _/A . 

dt dT 

/(r(t)) = /eg,A(r(0))exp [-/q AdT 


(23) 


= /eg,A(r(0))exp[(A)tt], 

where (.)* denotes the time-average. The term A signihes the phase-space compression factor, 
which denotes the average rate at which the phase-space collapses onto a fractal dimension 


smaller than the ostensible dimension 


29| . A can be related to several important dynamical 
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variables like Lyapunov exponem 


and entropy production (S) 


23 


25|: 


22| . and thermodynamic variables like heat flow {Q) 301 


S = -^ = fe(A), 


(24) 


The phase-space compresses (or expands) due to the heat flows from the individual ther- 
mostatted regions (the intermediate regions do not contribute towards phase-space compres¬ 
sion owing to Hamilton’s evolution equation), and may be split up into two parts: 


A = Ah + Ac = 


Qh Qc 
ksTn ksTc 


(25) 


Utilizing equation ((23]), we may write the nonequilibrium distribution post time t as: 


/ne,(r(t)) 


/eg,A(r(0))exp 



/eg,yi(r(0)) exp [((Ah)* + (Ac)*) t] 


(26) 


Equation fl26ll represents the general nature of a nonequilibrium distribution function, and 
therefore, represents a steady-state distribution function as well. Another important con¬ 
clusion from equation fl26l) is: 


Ue,{T{t))dT{t) = /,,,^(r(0))dr(0) 


dT{0) 

dr{t) 


= exp 


- lo A*' 


(27) 


The generalized work function, AfUc, introduced in the main text is the same as the one in 


Williams et. al. 


26j, and relates two microstates that are not necessarily in equilibrium: 


exp [AWc] 


/i(r(o))dr(o)z,, 
Mm)dT{t)z,, • 


(28) 


For deriving equation (7) oi the manuscript, we look at the transition between the equilib¬ 
rium state A, and the nonequilibrium steady-state obtained after introducing the thermal 
gradient (for a time t). During this transition, A does not change, and as a result, the 
partition functions may be omitted. Here, /i denotes the canonical distribution function 
shown in equation (|22D . and /2 denotes the canonical distribution function associated with 
the nonequilibrium microstate T{t). Appropriate substitution results in: 


exp [AHc] 


exp[-/3oE(r(0))]dr(0) 
exp [-/3oE(r(f))] dT{t) ■ 


(29) 
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The ratio of the differential volume terms are related to the phase-space compression factor: 


dT{0) 

dr{t) 

Substituting equation fl30|) in equation 


= exp 


Adt' 


(30) 


, we get: 


tattat exp [—/3o-E(r(0))] 
exp [AITg] = __ r o ^ ^Xp 


Adt' 


(31) 


exp|-/lo-B(r(i))] L jj 

Employing the first law of thermodynamics, equations flM|) and fl25l) . and recognizing that 
no external work is performed during the transition from A —)■ S'S'i, we can write: 


E = Q 


Substituting equation 


dt' 


E{r{t)) = E(m) + /„■ + Qc 

= £(r(0)) + kB J‘ ITh{a„) + rc(Ao)] dt'. 

into equation ([3T]), we get: 


(32) 


exp [/AWg] = exp [-/3o (^(r(0)) - E{T{t)))] exp - Adt' 


exp [AITg] = exp 

<3o(/o‘(Qh + Qc)*') 

exp 

exp [AITg] = exp 

^0 (/o Qtj + Qc) dt'^ 

exp 


~ Jo 

_ L fM _i_ 

Jo I 


dt' 


(33) 


exp [AfTc] = exp 
exp [AITg] = exp 


/^o/c 


t I Th — Tq ^ ^ Tc — Tq 

0 


-Qh + 


Th ^ Tc 

^ - ^Qc,t 


Qc ) dt' 


Th ^ ' Tc 

Since the relation fl5^ holds true for a generalized nonequilibrium state, it must hold true 
for the nonequilibrium steady-state as well. Therefore, we write: 

= A (^Qh., - A^Qc. 

\ -lh J-c 

which is the same as the equation (7) of the manuscript. 


t , 


B. Section-II 


We will use the generalized work function to (i) derive the fluctuation theorem for heat 


fow 


3l| . and (ii) calculate the probability of violation of Fourier’s law in thermal conduction 
30j |. Rewriting fl33|) in terms of the phase-space compression factors, we get: 

rt 


AITg = f^o ( I [kBATnAn — ksATcAc] dt' ) , 


(34) 
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We now take the special case that Tq is the average of Th and Tc 

^»+Z£;Ar„ = Arc = hL-^c 


30 


m 


To = 

Snbstitnting this in equation 
1 


2 

, we get: 


AWg = 


Th + Tc 


[{Th — Tc)Ah — {Th — Tc)Ac] dt' — 


'-H 


’-C 


Th + Tc 


[Ah — Ac] dt'. 


(35) 


For a Nose-Hoover therniostatted system in a d-dimensional phase-space, having Nh par¬ 
ticles under the influence of the thermostats, the phase-space compression may be written 
as: 

Ah = —dNTTjH'i Ac = —dNTTjc, (36) 

where rj represents the Nose-Hoover reservoir variable. Equation (135|) may now be written 
as: 


A ^ Ar f 

AWc; - dNx 

.LH + .Lc Jo 


[Vc - Vh] dt'. 


(37) 


We have introduced the superscript F to denote the time-forward motion. Equation (l37|) is 
the same as the equation IJ derived by Evans et. al. [30||. It is evident that AWc depends 
upon the initial microstate from which the trajectory initiates. Therefore, AW^ in equation 
(137|) may be written as AfF^ (r(0)). To obtain the fluctuation theorem for heat flow 3l(|, we 
will look at the time-reversed dynamics. In the time-reversed dynamics, the system begins at 
the nonequilibrium microstate r*(0) which is the same microstate as r(t) but with reversed 
momenta. The system reaches in time t the microstate r*{t) which is the same microstate 
as r(0) but with reversed momenta again. Therefore, the energy functions may be related 
as: 


E|r*(i)l = -B|r*( 0 )]+ % 

E ir(o)| = E |r(()| + /' \q-h + Qh 


dt 

dt 


(38) 


Qh + Qc 


-- Q*H + Qh 

While writing the last equality, we have used the relation fl32l) . In simple terms, the equation 
(138|) says that the heat flow from the thermostats in the time reversed dynamics is exactly 
equal and opposite to the one in the time-forward dynamics. The generalized work function, 
therefore, during the time-reversed transition becomes: 

ATc 


A»'«(r*(o)) = 


T, 


C 


Qc.t — /^O 


^QH,t - ^Qc,t 
Th Tc J (39) 


Aw|(r*(o)) = -Aw|(r(o)), 
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Since there is no change of A during A —>■ SSi, the time-reversed trajectory represents the 


conjugate trajectory moving forward in time. In the same terminology as Evans et. al. 30| . 
one may therefore, view AIEQ(r*(0)) synonymously with AlE(|’(r*(0)). Now, we relate the 
probability of observing a trajectory to its conjugate trajectory: 

p{AwS^-k) = /,,i(AH^^(r*(o)) + «:)/2(r*(o))dr*(o) 

= /ryAH^^(r-(o)) + fc)/2(r(())<ir(() 

= /rhAW"G(r’(0)) + k) exp [-AiF|(r(o))] /, (r(o))<ir(o) 

= y i (AH^'-(r(0)) - k) exp [-A»y|(r(0))] /, (r(o)) <ir(o) 

= exp [—k] P {AWq = k) 

Substituting AWq from equation (I57|) . we get: 

P [AW^ = k) 


P {AW^ = -k) 




J F / 


1-1 

SO 

1 

o 

. 

L 

1 + 

E^ E? 

A 

1_1 

= -fc) 

k J 


exp [k] 
exp [k] 



^ /o [Vc - Vh] dt 

= k) 

J F / 


Jo [vc - Vh] dt 

= -k) 
F J 


P{[r]c-r]H]F = k) 
P{\f]c-nH]F = -k) 


= exp 


= exp 


dN^ 


H 


-C; 


Th + Tc 
dNr ^J^ ~ kt 


Th + T, 


which is exactly what is derived by Evans et. al. in equation (15) of 301. 


c 


(40) 


(41) 


Section-Ill 


In this section we derive equations (8) - (12) of the manuscript. Proceeding analogously 
like in the previous sections of the supplementary material, but now realizing that during 
the transition period r, the work done also features in the first law equation fl32l) . we can 
write: 


rt+T 

E{T{t + r)) = E(r(0)) + W + kB [Th{Ah) + Tc{Ac)] dt'. 

Jo 


The generalized dimensionless work function now becomes: 

exp [-/3oE(r(0))] 

exp [AW^^a^ss 2 ] = 




dr(o) 


exp [-^oE{T{t + r))] 


■'Xb 


dT{t + 


T 


X 


^Xa 




(42) 


(43) 
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The partition function at t + r is Z\g because its associated equilibrium state is B. The 
equation gets simplihed into: 


exp = exp [-/3o (^(r(0)) - E{T{t + r)))] x exp - £ Adt' 


exp [j^oW] X exp 
exp [(3qW] X exp 
exp [(3qW] X exp 


/< 


t+r 


T.-'h-ii-"-'’ 


^ si*’ |Ar„A„ - ATcAc] dt' 

-1-0 


f^O Jr 


t+r 


^Qh - ^Qc 

j- i 


H 


■C 


dt' 


dt' 


(44) 


Therefore, the generalized work function becomes: 


exp [AVP|,^^ssJ = m + A 


B'H,t+T 7 ^ — ^C,t+T 

-L I 


T, 


H 


C 


(45) 


which is same as equation (8) of the manuscript. Subtracting equation fl33|) from equation 
gives the following: 




ffoW + A 
^oW + P„ 

feq,A (r55j^ 


^ {QH,t+r - QH,t) — 7f;^{Qc,t+T - Qc,t) 

1 ( 


Th 

ATh ATc 

- y —VC, 

0)) dr55^(0)ZA^ 


c 


feq,B ('?■)) dTsS 2 {T')ZxB 


(46) 

By looking at the dehnition of generalized work function, it is evident that the last equality 
of equation (1461) gives the generalized work function during SSi SS 2 i.e. 


AW^C ,- AhT|,A^55, = m + /5o 

^ /eg,A (r5gi(0)) drg5^(0)ZA^ 
feq,B i^SSiiT)) dV sS 2 {'^)Z\b 

/eg,A(r 55 i( 0 ))dr 55 ^( 0 ) represents the probability of the nonequilibrium microstate rss.(o) 
in the associated equilibrium state A. Likewise, feq,BiX 882 ( 7 ))ss 2 i.^) is the probability in 
the associated equilibrium state B. Equation fl47|) is the same as the equations (9) and (11) 
of the manuscript. In a similar manner, by looking at the reverse transition, one can derive 
the equation (10) of the manuscript. We will now drop all subscripts except G. One can use 
ergodic consistency - every microstate in steady states SSi and SS 2 can be obtained from 


- ^Qc,d 


'-H 


T, 


c 
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equilibrium states A and B, to derive equation (12) of the manuscript: 


P (AW^ = -k) = /r ^ (AWf^ + k) feq,B {'^sS2{t)) dTss2{'r) (by using Liouville’s equation) 


= /r 5 {AW§ + k) exp [-AW^] Ua (rs5i(0)) dVssA^) 
= /r <5 {AWA - k) exp [-AWA] feq,A (r5Si(0)) dTssA^) 
= exp [-k] {AW^ = k) 


Zxs 

Z\a 

Zxb 


(48) 

Ideally, one should be using the nonequilibrium distributions at time 0 and r, and not the 
equilibrium distributions while deriving the previous expression. However, because of a lack 
of such nonequilibrium distributions, we are limited to using equilibrium distribution func¬ 
tions. As a consequence, the contributions of phase-space compressions, which seep into the 
dynamics when nonequilibrium conditions are imposed, cannot be accounted. Therefore, our 
method works only for large r. Taking r large enough, while hxing the time required to reach 
the steady state, ensures that the contributions arising from the phase-space compressions 
become negligible. 
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